Calculation of the one-particle and two-particle condensates in He-II at 

r = 



^ ' Maksim Tomchenko 
O 

o 

. Bogolyubov Institute for Theoretical Physics, 

^ I Metrologichna St. 14-b, 03143 Kiev, Ukraine 

E-maih mtomchenko@bitp.kiev.ua 

in 



5-H 



We analyze the microstructure of He-II in the framework of the method of collective variables 
(CV), which was proposed by Bogolyubov and Zubarev and was developed later by Yukhnovskii and 
QJ Vakarchuk. The logarithm of the ground-state wave function of He-II, ln\l/o, is calculated in the 
approximation of "two sums", i.e., as a Jastrow function and first (three-particle) correction. In the 
P CV method equations for \l/o are deduced from the N -particle Schrddinger equation. We also take into 
', account the connection between the structure factor and \l/o, which allows one to obtain \l/o from the 
^ structure factor of He-II, not from a model potential of interaction between He-II atoms. It should 
I ■ be emphasized that the model does not have any free parameters or functions. The amount of one- 
: particle (Ni) and two-particle (N2) condensates is calculated for the ground state of He-H: we find 
O , Ni ^ 0.27N and N2 ~ 0.53A^ in the Jastrow approximation for \l/o, and, taking into account the 
i three-particle correction to ln\l/o, we obtain Ni ^ 0.06A^ (which agrees with the experiment) and 

^ ■ N2 ~ O.lGiV. In the approximation of "two sums", we also find that the higher s-particle condensates 
>- (s >^) are absent in He-II at T = 0. 
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The sum of one-particle condensate (IPC), two-particle condensate (2PC) and higher s-particle con- 
densates is usually referred to as the composite condensate. Knowing the structure of the composite 
condensate in He-II is, without doubt, of great importance pQ-fH]. Except for a purely cognitive 
interest, it also has a "practical" side, namely, in the field-theoretic approaches to the modelling of 
He-II microstructure [U SI El [11] , the quasiparticle spectrum of He-II is explicitly expressed through 
the amount of the IPC and 2PC, and the dependence on higher condensates is also not excluded. Of 
interest is the question whether all of He-II atoms belong to the composite condensate at T = 0, as was 
suggested in [ZllHlITT]. The superfluidity of He-II by itself is probably [3l[T2] caused by an off-diagonal 
long-range order (ODLRO). 

Our work is devoted to a calculation of the amount of IPC and 2PC in He-II at T = 0. To describe 
the microstructure of He-II, we use the method of collective variables (CV), which was first proposed 
by Bogolyubov and Zubarev and was later developed in the works by Yukhnovskii and Vakarchuk 
[T1]-[T8]. First, we obtain the ground-state wave function of He-II, \E'o, and then we calculate the 
amount of the condensates using the formula of [T71 [19] for the s-particle density matrices Fg. The 
model does not contain any free parameters or functions: \E'o is obtained as an eigenfunction of the A^- 
particle Schrodinger equation; we also take into account the connection between and the structure 
factor S{k) of He-II 
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As far as we know, the higher s-particle condensates {s > 3) were not calculated previously. The 
amount of IPC in He-II at T = was found in many works 0, [H [ISl [2n]^[21l [H] (though only in 
[T8l [23t [23] it was done without free parameters) , and the theory agrees with the experiment on 
the whole. The amount of 2PC in He-II is still unknown. The IPC was measured in many works, as 
a number of atoms with momentum equal to zero, but how to measure the 2PC is not yet clear (this 
issue is beyond of the scope of the article). There are only several theoretical estimates of the 2PC 
[HI [8], all of which use free parameters. The results of [6] are discussed below in Sec. 4. In [8], a final 
result was not presented, but some equations were obtained, from which it follows that the proportion 
between IPC and 2PC can vary depending on the form of the potential of the interaction between 
He^ atoms. We calculate the amounts of IPC and 2PC without free parameters, and in more exact 
approximation as compared with |TH HH [6] . 

In and \l8l, the amount of IPC was found in the approximation of "one sum" (IS) with the 
result A^i ^ 0.08A^ and Ni ^ 0.04A^, respectively. Below, we obtain IPC (A^i) and 2PC (iVa) for He-II 
at T = in more exact approximation of "two sums" (2S), and we find A'^i ~ 0.06A^ and N2 ~ 0.16A^. 
Our formula for 2PC refines the Ristig's formula [6] previously obtained by another method in the 
IS-approximation for ln\E'o- 

In [19j we have found also that (i) in a weakly interacting Bose gas, all atoms belong to IPC or 
2PC at T = 0, and (ii) the higher s-particle condensates (s > 3) are absent in He-II at T = (which 
was shown in the 2S-approximation). 

2 The ground-state wave function of He-II 

There exist several methods for calculation of the ground-state wave function of He-II (see also the 
review |25]): the variational method |26j, the "Green's function Monte Carlo" method [27j and it's 
development, "shadow wave function" (SWF) method [28], the Path Integral Monte Carlo simulations 
[22], the diffusion Monte Carlo (MC) simulations [231 the "hypernetted chain" (HNC) method 
|29t [50], (all of these are indirect methods for solving the iV-particle Schrodinger equation), Feenberg's 
approach [21] (solving the Schrodinger equation in the r-space), and the CV-method [I3]^[in], [321 E3] 
(solving the Schrodinger equation in the k-space). 

The quantum-mechanical models are developed actively and a noticeable progress is already 
achieved [HI [321 [231 [25l [281 [301 [33] • Modern variational methods [28] (SWF), [29ii30j (HNC) reach an 
accuracy of the order better then 0.1 ii', but such models, unfortunately, use several free parameters. 

In the CV-method, the equations are deduced from the first principles, namely, from the exact 
A^-particle Schrodinger equation, and a solution of these equations can be found numerically without 
using any fitting parameters or functions (which is important) by taking into account the connection 
between \I/o and S{k). 

All existing models of He-II of which we are aware (except the mentioned CV-approach) use several 
fitting parameters, at least, in the effective interaction potential (even the MC simulations); in this 
case, of course, it is not so difficult to obtain two "points": the ground-state energy and IPC. In 
our opinion, the approach without fitting parameters is preferable because, in this case, nothing is 
introduced in the model "by hand" . In our paper, we use the CV-method since this method does not 
have any fitting parameters. 

According to [15], the ground-state wave function of He-II has the form 

^0 = -y^, So = Y^ — - — '^(ki + • • • + kj)aj(ki, . . . , kj)pk, . . . pk., (1) 

VV j>2 J- ki,...,kj5^0 

where 6 is the Kronecker delta, is the total number of atoms in helium, Q is a normalization 
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constant, and pk are the collective variables 



1 ^ 



Pk = ^Ee-^''■^ k^o. (2) 



The term with a„ in ([1]) corresponds to n-particle correlations in r-space. The zero-order approximation 
for \l/o, which is also the IS-approximation, is 

1 ^ x-^ CLoik) , , 

In $0 = In Q + 2l ^-^PkP-k, an>3 = 0. (3) 
The wave function \l/o ([3]) can be present in the well-known Jastrow form 

^o = ^I[e''^-^"'\ (4) 

where 

= l5:^e^'^r = /rfk^e^'^^ (5) 

^ ' ^ k 2 Ar(27r)3 7 2' ^ ^ 

and V is the volume of the system. 

The capabilities of computers force us to restrict ourselves to the 2S-approximation for \l/o; in 
this case, the sums with 02 and 03 are taken into account in ([1]), but a„>4 = (for the "n-sum" 
approximation [TTJIIB] all sums up to the n-th sum over k are included in all expansions). Substitution 
of \E'o dl]) into the iV-particle Schrodinger equation allows one to obtain a chain of equations for 
a„. In the 2S-approximation, we have 

(u u \ _ 2a2(fci)a2(fc2)kik2 + 2a2(fci)a2(fc3)kik3 + 2a2{k2)a2{h)k2k3 
""'^ kl[l - 2a2ik,)] + kill - 2a2ik2)] + kill - 2a2ik,)] " 

Throughout in the paper, we assume k3 = — ki — k2 7^ 0. The s-particle density matrices Fg for the 
ground-state of Bose liquid were found in [17] for the approximation of IS, and in [TTJ [19] for the 
approximation of 2S . Moreover, in [I7l[T9], an equation connecting \l/o with the structure factor S{k) 
was obtained: 

2a(k)-l-^- ^^(^) - ^ (7) 
Sik) Sik)[l-2a2m Sik)' 



where 



y (k) = — a2(g)a2(k + g) + a3(k, g) [1 + a3(k, g)] 
N^^, [l-2a2(g)][l-2a2(k + g)] ' 



iV^^Jl-2a2(g)][l-2a2(k + g)]' 
and also the formula for the amount of IPC (A^^i) was found: 

In (iVi/iV) = IiA + I2A + I2B + he, (10) 



8iV2 i ,11 1 _ 2a,(k,) j 1 - 2a2(fci) ' ^ ^ 



'2_B 



"2/V2 X! I n 1 _ 

ki,k27^0 \j=l ^ 

2a2(fci) + a3(ki,k2) 
1 - 2a2{ki) 



2a2{kj 



a3(ki,k2) * 

2a2(A;i)[l- 02(^:2)] J 
1 



(13) 



(14) 



Note that equations similar to (j6]), (171) were obtained by other methods by Campbell and Krotschek 
[29] . Our system of Eqs. ([6])-([9]) is a little more exact, because in (|6]) we have with 02 from ([7]) 
(which includes correction 03 7^ 0), but in [29j 03 was taken with 02 in zeroth approximation (03 = 
in (HD). 

The following asymptotics as A; ^ is true |19j : 



2a2(A; ^ 0) 



1 + ^2(0) 
Sik) ' 



21,2 



E{k ^0) = ck 



S{k^O) 



h'k 



2mS{k) 
Tik 



:i + S2(0)) 



2mc 



(1 + ^2(0)), 



where E{k) is the quasiparticle spectrum of Bose liquid (see the equations for E{k) in [16 
and 

02 (g) 



(15) 

(16) 

(17) 
[33]), 



S2(0) = -E 



n 4 



l-2a2(g) 



> 0. 



In the IS-approximation ([3]), for which 

2a2(A;) 



S{k) 



fl?i>3 — 0, 



(19) 



for He-II we obtain ^2(0) = 0.33 (from ( fT8|T9l) ). and, in the 2S-approximation, a numerical solution 
of (jnD-dS]) (see Sec. 3) gives ^2(0) = 0.66. 



3 One-particle condensate in He-II at T = 

In the IS-approximation, the amount of IPC can be simply found from ( ITOl) . ( ITTi) . and ( |T9i) . For the 
structure factor S{k), we use the smoothed experimental data from [36j, which is, perhaps, the most 
exact. We extrapolate this data to T = according to 



S{k, T = 0) = Sik, T) tanh (20) 

At T = 0, we should have the asymptotic S{k = 0) = |38]; in fl20|) . we take into account that 
E{k 0) = ck. From ([lO]), ([II]), ([ED, and ([20]), we numerically find rii = ^ ■ 100% = 27.2%. 

To calculate A^^i in the 2S-approximation, we need to know, according to f|IO]) - f|I^ . the functions 
02 (fc) and a3(ki, k2). The function 03 is defined in and to obtain 02 (A;) we should solve numerically 
the integral equation (I7j) taking into account ([6]), ([H]), and ([9]). We everywhere replace sums by integrals 
according to the rule [IS], [U] 
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Equation ([7]) cannot be solved by iteration; we have succeeded in solving it numerically by the Newton's 
method, in which way we found a2{k). From (jH]), (niI|) - f|T^ . and fl^Tl) . we obtain rii = 6.1%, which 
agrees with the experiment: nl^^ ^ 6 — 12% for T = p9j-[41j. We estimate the numerical error for 
ni to be — ~ 0.1. The properties of ^'o and the possibility to obtain the ground-state energy £"0 are 
discussed in more detail in [32] . 



4 Two-particle condensate 

It can be shown in 2S-approximation that the s-particle density matrices Fg of helium-II for T = in 
three dimensions display ODLRO: 



lim Fs(ri, . . .,rs\r[, . . . ,r'J|ij = F^(oo) = [Fi(cx))]'* = const > 0, 

|ri-r'-|^oo 



(22) 



B: |rj — Fjl, |r^ — r^l are fixed for any i,j = l,...,s, i ^ j. 

From the general principles, we write the probability M/ki,k2 of finding the momenta ki and k2 in two 
arbitrary atoms, for T = 0: 



= ^ / ^2(ri, r^K v'^y''^^''^-<^^'''^^'^-''^'^dv^dv2dv[dv'^ 



(23) 



where F2 is the two-particle density matrix. 

To obtain the amount of 2PC, we should know H^k.-k for T = 0. Let us consider first the IS- 
approximation for and F2. In this approximation, is defined by ([3]), and F, = F]^) was found 
in [HI: 



FW(ri,...,r,|r;,...,r- 



exp I 

s a2{k) 
iVl - 2a2{k) 



where 



h{k)=a2{k) + f2{k), f2{k) 



+ f2{k)^^.C 



am 
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l-2a2(A;)' 



-ikr' 

e J . 



(24) 

(25) 
(26) 



It is convenient to represent ^2^^ in the form 



InF^'^ = lnF2(oo) + ^i(ri - r2) + v^i(r; - r'2) + 
+ ip2{r[ - ri) + V92(r2 - r2) + ^2(^2 - ri) + Lp2{r[ - r2] 



The functions have the property (pi{r 00) ^ ^ —>■ 0. From fl23|) . we obtain 



k,-k 



^2(00) 



I rfric/r2c/r;dr;,e^''("^-"'i)-*''('"2-'"2)exp [y?i(ri - ra) + 



(27) 
(28) 

(29) 



+ ipi{r[ - r'a) + (p2{r[ - ri) + v?2(r2 - r2) + ip2{r'2 - ri) + (p2{r[ - r2] 



To calculate VFk,-k (EHD, we expand the exponent in 
can be shown, that only integrals of terms of the form 



with the sum of all yjj, in a power series. It 



tp^l (ri - (r'l - r^) and {v[ - n)'^^^ (r'2 - ra) 



(30) 



[Ij = 1, 2, 3, . . . are all possible natural numbers) are significant for the value of W^k,-k- The integrals 
of other sets of ifi are times smaller or equal to zero. Taking all this into account, we obtain 

^2(00) 



k,-k 



$,(fc) = -J e'^'(^)e*'^^rfR, F2(oo) 
Following the notation of [6], we rewrite PVk-k (Ell) in the form 



where 



iV^l^k.-k = iVkiV_k + XkX-k, A; ^ 0, 



iVk = Fi(oo)<l>2(A;), Xk = Fi(oo)<l>i(fc). 



Here, A^k is the amount of atoms with momentum k. For A; = 0, we have (f2{'r) — >■ as r 
$2(0) = N] therefore, iVk=o = NFi{oo), and we obtain the one-particle condensate Ni 



(31) 
(32) 

(33) 

(34) 
00, and 



k=0- 



In 



^2(00) = F^^(oo), and Fj(oo) correspond to the IS-approximation: ^ = Fi(oo) = e^^^ 



dSH) and 
see (HUD, (dH). 

As it can be seen from (1331) . W^k.-k consists of two terms: the first one, j^, is simply the product 
of relative numbers of atoms with momenta k and — k, this term does not describe correlations. The 

I |2 

second term, describes correlations in k-space in the pairs (k, — k). It is naturally to relate the 

two-particle condensate precisely with the correlation term 

The 2PC was already calculated by Ristig [6] using a different method. In Ristig's works [6], the 
2PC was determined by the quantity 



k^O 



k^O 



(35) 



The expression P2 can be interpreted as the mean degree of correlation of the pairs of atoms with 
momenta ki + k2 = 0, k, 7^ 0. We believe that it is more reasonable to define the number of atoms in 
2PC as follows: 

A^2 = E ^kCk, (36) 

k^O 



where Ck is the "correlation factor'' 



Ck 



IXkl 


for 


IXkl 


< ATk 


1 


for 


IXkl 


> ATk 



(37) 



Using 02 (fc) in ( fTOll . from (l25ll -( IMl) we obtain the functions Xk and A^k (see Fig. 1); then from ( |35l) - (l37j) 
we find A"2 ~ 0.53A^ and P2 ~ 0.31. In [6], using the Jastrow approximation (jl]) for \E'o [as in our work 
([3])], it was found P2 ~ 0.09, but, for S'i(r) in Ristig applied P the popular McMillan's form 



ro = 2.963A. (38) 
This form is simple and is convenient for calculation but, at the same time, it is a very crude variational 




Fig. 1. The functions A^k and Xk, where k is in 
the units of fc^ = 2Tr/d = 1.756A"^. Open cir- 
cles mark iVk in the IS-approximation; triangles mark 
A^k in the 2S-approximation; crosses mark Xk in 
the IS-approximation; squares mark Xk in the 2S- 
approximation. 



Fig. 2. The function 5i(r/a) ©, where a = 2. 641 is the 
diameter [42] of a He''' atom. Open circles correspond 
to the IS-approximation for a2(fc); crosses to the 2S- 
approximation for 02 (fc); triangles to the "three sums" 
approximation for 02 (A;) in the approach of |35| with a 
model elliptical potential with U{r = 0) = 60K; contin- 
uous line shows the McMillan's form ([38| . 



approximation. Using the CV- method, we find Si{r) (see Fig. 2) significantly more accurately from 
([7]) and ([5]) without any free parameters. Equation (I7j) was deduced [HI [19] from the equation 

S{k) = l + yj dv{g{r) - l)e-^''^ (39) 

connecting S{k) with the pair distribution function g{r). Equation ([7]) was derived in the approx- 
imation neglecting a„>4 in ([1]) and neglecting multiple scattering in fl39|) |33]. As shown in Fig. 2, 
our result 5*1 (r) differs greatly from fl38l) : therefore, our result for P2 differs from that obtained in [6] 
although the formulas for Xk [see (1451) below] and P2 are the same in our work and in [6]. Note that 
Si{r 00) ~ ^ since a^^k ^ 0) ~ -i, see ([7]). 

It should also be noted that we found the function Siijr) not quite presicely at r < (ic ~ 2 A (at such 
r, two atoms penetrate each other). This is evident from the fact that the function g{r) recovered from 
Eq. (139|) oscillates irregularly at r < 1.7 A (see Fig. 3) and is different from that expected from physical 
considerations: g{r) —> at r < cic, and g{r) > 0. Inaccuracy of the determination of g{r) and Si{r) 
at r < dc is caused by large relative error of measurement of the value S{k) — 1 at A; > 27r/(ic ~ 3 A~^, 
probably by neglecting multiple scattering in ( 139|1 [13] . and by truncation of expansion ([T]) [for 5*1 (r)]. 

Oscillations of g{r) at small r arouse mainly [36] from the inaccuracy of measurement of S{k) at 
large k, k > 27i/dc- Although the values of S{k) are very close to unity at such k, S{k) ~ 1, the 
values of g{r) are sensitive to the small oscillations of S{k) in the neighborhood of the unity. At the 
same time, such a small indefiniteness of S{k) at large k infiuences insignificantly the amounts of the 
condensates and quasiparticle spectrum [32]. According to our numerical analysis, the amount of the 
condensates are sensitive, first of all, to the value of S{k) at the small and middle k, k <2iT/d ^ 1.7 A~^ 
(where d = 3.58A is the mean distance between helium atoms), which corresponds to r > d {g{r) is 
well defined at such r). In fact, the condensates are "spread out" in the whole system, so the values 
of g{r) at large r, r > d, are important for the estimates of the amount of condensates. Moreover, the 
oscillations of g{r) at small r may be related to the problem of realistic description of the structure of 
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0,0 0,5 1,0 1,5 2,0 2,5 3,0 
r/d 



Fig. 3. The pair distribution function g{r) of He-II at T = 0, from p9l 



helium atoms. 

Correlations in the pairs (0,0) are absent |19j; therefore, only atoms with A; > belong to the 
two-particle condensate (k, — k). 

In the more exact 2S-approximation, the two-particle condensate is calculated similarly (see for 
more details). The function I^ic-k is again defined by (13T1) - (!MI) with ipi given by (128|) . where Fj(oo) 
now correspond to the 2S-approximation [see ([I0])-([T1D for ^ = Fi(oo)], and we derive more exact 
formulas for fi{k) and f2{k)'- 

(40) 
(41) 

(42) 
(43) 



N 

fiik) = a2{k) + f2ik) + 6i{k), 

1 ^ ct3(k,q) 

1 _ 

q^O 



Nl-2a2{k^)^j^J}^l-2a2{k,y 



ka^O \j=l ^ 



2a2{kj 

a2{ki) + asiki, kg) (1 - a2iki)) 



a3(ki,k2) 



+ 02(^1) (1 - 202(^2)) + a2(fc2)a2(ki + k2) 



(44) 



1 - 2a2(A;i) 

The function Xk in the 2S-approximation can be presented in the form 



(45) 



where Fi{r) = ^e*^^ (this is the one-particle density matrix in the 2 S- approximation, see Fig. 4), 
S'i(r) is defined according to ([5]), and 



S*2ir) 



N 



,jkr 



(46) 



S 




Fig. 4. The one-particle density matrix Fi{r) for He-II. 
Circles correspond to IS-approximation, crosses corre- 
spond to 2S-approximation. The dotted line shows the 
amount of one-particle condensate, ^ = Fi{oo), for the 
approximations of IS and 2S. 



Fig. 5. The fmiction x{r)/nc (gT]), where d = 3.5781 is 
the mean distance between He-II atoms, and ric = Ni/V. 
Circles and crosses correspond to the same approxima- 
tions as in Fig. 4; the continuous line shows the approxi- 
mate form of the functions g{r) and x('')/^c for a weakly 
interacting Bose gas [19]; the dotted line indicates the 
level x('')/^c = 1- 



For the IS-approximation, we have 6i{k) = 0, 5*2 (r) = 0, and Fi{r) = Fi (r), so that ( l45l) coincides 
with the Ristig's formula for Xk derived earher [6] (to an accuracy of the factor N/V appearing as a 
resuh of different normahzation of Fi, we also choose the opposite sign for Xk)- Formula fHS]) was not 
obtained previously. 

Using a3(ki,k2) given by ([6]), and a2{k) as a solution of ([7j), we find the two-particle condensate 
according to ([MM), (gOD-dlH]): N2 ^ 0.16A^, P2 ^ 0.025. The obtained functions Xk and A^'k are 
shown in Fig. 1. 

Note that, in the 2 S- approximation, we have, taking into account (fT5|) . iVk^o = "^^^^ = 
+ ^2(0)), Xk^o = —Nk^Q. Such an asymptotics for A^k refines the one found earlier in 

[Sl@S],iVk^o = fi|. 

In Fig. 5, we show the function 

X{r) = J^J ^ke-^^^Xk = ^e^^^«+^5«Fi(r) = ^e^^W. (47) 

The function g{r) describes the dependence of the correlations in the pairs of two atoms with arbitrary 
momenta, (ki, k2), on the size of the pair r; and x{^) describes the same for a narrower class of pairs 
(k, — k). For infinite separation between the atoms in the pair, correlations disappear; therefore, g{oo) 
and x(oo) take the values of g and x for uncorrelated pairs [for a totally chaotic distribution of atoms 
in the r-space, we would have g{r) = g{oo) for all r, and similarly for x(^)]- 

Comparing the functions g{r) (Fig. 3) and x{^) (Fig- 5), we can see that these functions have similar 
form, as it could be expected. Atoms interact strongly at a small distance, so the correlations in the 
pares are significant at r < 2d, where the values of g{r) and x('^) are substantially different from those 
at infinity. 

The authors of [2T1 [TO] believe to be warrant the investigations of the possible pairing of He-II 
atoms in the r-space. However, in our opinion, the properties of x{^) and the results of [21J do not 
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give grounds to talk about the bound states in the r-space; these properties are indicative only of 
correlations in pairs. The pairing energy Ep ^ —0.05 K was found in [2T], but this energy is too 
small for the appearance of real bound pairs because Ep is two orders of magnitude smaller than the 
ground-state energy Eq. 

At present, we are not aware of any reliable evidence of the existence of bound pairs of atoms in 
He-II. An interesting arguments in favour of the possibility of pairing were presented in [T0|. Our 
analysis cannot exclude the possibility of existence of narrow bound pairs, of atomic size. Electrons 
should be collectivized in narrow pairs; therefore, such pairs have to be considered as a new kind of 
"atoms"; in this case, we need to know \I'o for a mixture of two Bose liquids. However, the closeness 
of the experimental value of circulation [IB] to the value of one quantum, k = h/m, indicates that the 
number of bound pairs in He-II is very small or equal to zero. 

5 The higher condensates 

We also investigated correlations (in the k-space) in groups of s atoms {s > 3) under the condition 
ki + . . . -|- ks = for momenta kj of atoms. We were interested in the "nontrivial" part of the s- 
particle condensate, which cannot be reduced to the "lower" condensates; so we put ki -|- . . . -|- k; 7^ 
for all / = 1, . . . , s — 1. There is a wide-spread, but not proved, opinion according to which the higher 
condensates should be present in He-II. We considered the problem in the 2S-approximation, for which 
the three-particle term with 03 is taken into account in ln\E'o dlD? ^^^d we expected that at least three- 
particle condensate would turn out to be nonzero. However, we found [inj that all higher s-particle 
condensates (s > 3) are absent from He-II at T = 0. An analysis is cumbersome and similar to that 
for 2PC, see [TS] for more detail. 

Our results, which consist in the absence of all higher s-particle condensates (s > 3) and in the 
small value of two-particle condensate, should be important for the field-theoretic model of He-II: 
they can significantly simplify such model and, perhaps, open a possibility of constructing an "ideal" 
micro-model of He-II. 

6 Conclusions 

The main results of our paper were already outlined in the introduction. Summarizing, we notice that 
the CV-method allows one to calculate the quasiparticle spectrum of He-II [El [32l [35l [33] and also, with 
lower accuracy, condensates (this work), without any fitting parameters or functions. The equations of 
the model are deduced from the exact microscopic equations, and the results approximately agree with 
the experiment. Thus, the model under consideration, undoubtedly, gives an approximate description 
of the microstructure of He-II. 

However, the calculated amount of one- and two-particle condensates appreciably depends on the 
number of the corrections an taken into account in ln\I'o. The reason is that the quantity ln\l/o and 
the condensates Ni are expanded in series in a parameter which is not very small (this is the function 
q{k) = a2{k)k/kd, kd = 27i/d, where d is the mean distance between atoms; the mean value of q{k) at 
k < kd is not very small, being around —1/2, see Fig. 1 in [19]) and, for condensates, the series stand 
in the exponent. These difficulties arise in all approaches where condensates are calculated from "^q, 
unless the free parameters are chosen so that the corrections "turn out" to be small. 

If, in some quantum mechanical approach, the amount of the condensate changes insignificantly 
by taking into account the next corrections to In \E'o and the condensates, this is, probably, mainly a 
result of a good choose of fitting parameters. Really, in our model we find a solution from the exact 
(but truncated) microscopic equations, without free parameters whatsoever, and we clearly see that the 
corrections to In \E'o and especially to the values of the condensates are not small enough, unfortunately. 
Nevertheless, our results give the approximate estimates of the fraction of the condensates (the amount 



in 



of IPC agrees with the experiment) and, at the same time, the significant value of the corrections show 
that quantum- mechanical approaches (except, perhaps, the MC methods) do not offer methods for 
sufficiently precise calculation of the amount of the condensates. 

The Monte Carlo method [23l [2^ is one of the perspective approaches for the calculation of the 
IPC and 2PC but, unfortunately, we do not know the error of the numerical definition of IPC and 
2PC for the MC simulations. 

The absence of small parameter is a general problem of virtually all (except the MC simulations) 
known approaches to the description of the microstructure of He-II. Perhaps, in some approaches, 
in the future, one will succeed in calculating the IPC, 2PC and higher condensates more exactly, 
using expansions in small parameters only and without fitting parameters and unjustified postulates. 
Unfortunately, by now, such an "ideal" micromodel of He-II is not constructed. And it is not clear 
even, whether it is possible. 

The author is very grateful to Yurii V. Shtanov for valuable discussion. 
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